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Abstract. In this research announcement we outhne the methods used in 
our recent proof that the Leech lattice is the unique densest lattice in R-^*. 
Complete details will appear elsewhere, but here we illustrate our techniques 
by applying them to the case of lattice packings in R^, and we discuss the 
obstacles that arise in higher dimensions. 



1. Introduction 

Given any lattice A C M" (i.e., a discrete subgroup of rank n), one can form a 
sphere packing by centering congruent balls at the points of A, with radius as large 
as possible such that their interiors do not overlap. Which lattice maximizes the 
density of this packing, the fraction of covered by the balls? This question arises 
naturally in geometry, number theory, and information theory (sphere packings are 
error-correcting codes for a continuous channel, as opposed to a discrete channel). 

The densest lattice is known for n < 8. In each of these dimensions, it is a root 
lattice: Ai, A2, ^3, -D4, -D5, Eq, E-j, or E^. Furthermore each is unique up to 
scaling and isometries. The books jCS| and |M] are excellent sources of information 
on lattice packings. (See [Hal IHTl O Ilc72l IV| for the original papers.) 

What may be the most interesting lattices are the Eg, root lattice in and the 
Leech lattice in R^** (see |E] for a beautiful introduction). They are connected with 
many other branches of mathematics, and are undoubtedly the densest sphere pack- 
ings in their dimensions. However, nobody has been able to prove their optimality. 
The Leech lattice has not even been known to be the densest lattice in its dimension, 
let alone the densest sphere packing (it is possible that there are sphere packings 
denser than any lattice packing). By contrast, Blichfeldt ^I] proved in 1935 that 
E% is optimal among lattices, and in 1980 Vetcinkin proved its uniqueness. 

In this paper we announce the following theorem: 

Theorem 1.1. The Leech lattice is the unique densest lattice in M.'^^, up to scaling 
and isometries ofM.'^^. 

Our method also yields a new proof that Eg is the unique densest lattice in K®. 
The proofs for and M^'* can be found in |UK| . Here, we will apply the techniques 
to the case of the hexagonal lattice in R^. That result is of no research interest, 
because it is quite easy to prove using much simpler methods than ours. Our goal 
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is simply to illustrate our approach in a case involving no technical difficulties. At 
the end of the paper we will outline the additional obstacles that occur in and 



p24 



Our principal technical tool will be the Poisson summation formula. Let / : 
R" — > R be a Schwartz function (i.e., / is smooth and all its derivatives are rapidly 
decreasing). We define the Fourier transform / of / by 

m - / /(x)e^"<^'*> dx. 



If A C R" is a lattice, then the dual lattice A* is defined by 

A* ^{ye R" : (x, y) G Z for aU x e A}. 

(If wi , . . . , Vn is a basis of A, then the basis u* , . . . , v* of R" that is dual to wi , . . . , 
relative to the inner product forms a basis of A*.) The Poisson summation formula 



states that 



^^(^) = vol(R»/A) 



For a review of the standard proof, see Appendix 

The most effective way to make use of a nontrivial identity such as Poisson 
summation is to choose / to make the two sides as different as possible. We will 
apply it using a function / such that f{x) < for sufficiently large, but f{t) > 
for all t. This technique was used in |CE| to prove upper bounds on the sphere 
packing density (see also These bounds appear to be sharp in R^ and R^'', 

which would solve the full sphere packing problem in those dimensions, but the 
apparent sharpness has not been proved. It appears that the only dimensions in 
which it is sharp are 1, 2, 8, and 24, and the one case in which it has been proved 
is R^ 

Our approach to proving Theorem II. II is to combine the analytic methods from 
|CE| with geometric and combinatorial arguments that are specific to lattices. As 
mentioned above, we will present the details here for the two-dimensional case. 
We wish to show that the obvious hexagonal packing is the densest lattice packing 
of disks in M^. For simplicity we will choose the following normalization of the 
hexagonal lattice. Let A2 be the lattice with basis v, w, where 

/2^ 



^ 
and 

The advantage of this normalization is that A2 is an even integral lattice. In other 
words, all inner products of lattice points are integers, and furthermore every lattice 
point has even norm. Here, |wp — |wp — 2 and {v,w) = 1. 

Let A C R^ be any lattice that is at least as dense as A2. We scale A so 
that vol(R^/A) = 1 (we refer to vol(R^/A) as the covolume of A). Unfortunately, 
vol(R^/A2) = \/3, so the scalings of A and A2 are not compatible. That will cause 
us no problems, and in fact each scaling will prove convenient at a certain point in 
the proof. 

Our overall plan is as follows: 
(1) Use Poisson summation to study the short vectors in A. 
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(2) Apply this information to see that A is, up to scaling, a small perturbation 
of A2. 

(3) Prove that A2 is a strict local optimum for the packing density. 

(4) Conclude that either A is the same as A2 (up to scaling and isometrics) or 
it is strictly less dense, which would contradict our assumptions. 

Before beginning, it is worth noting that the density of A simply depends on 
its minimal vector length (i.e., the length of the shortest nonzero vectors in A). 
Because vol(M^/A) = 1, there is one sphere per unit volume in space, and the 
density of the packing equals the volume of such a sphere. More generally, the 
density of any lattice L C M" equals 

vo\{Br) 
vol(M"/L)' 

where r denotes the packing radius and Br denotes a ball of radius r. The packing 
radius is half of the minimal vector length, because that is the radius at which 
adjacent spheres are tangent. Thus, our goal becomes understanding the minimal 
vector length of A. 

Because A is at least as dense as A2, its minimal vectors have length at least 
(4/3)^/*. To see why, recall that vol(M^/A2) = \/3 and its minimal vectors have 
length a/2- Rcscaling A2 by a factor of 3~^/^ makes the covolume 1 and yields 
minimal vector length (4/3)^^''. 

2. Applications of Poisson summation 
Define / : ^ M by f{x) = p/(27r|a;p)e-'^l"'l', where 

Pfiu) = 20812 + 756u + llOTu^ - 216u^ 
Then one can calculate that f{t) = pj(27r|tp)e^^l*l , where 

pj{u) = 20812 + 5940w - 2781^2 + 216u^ = (43 + 24u)(-22 + 3uf. 

We have /(O) = /(O), f{x) < for \x\ > 1.084, and f{t) > for all t. (Terminating 
decimal expansions such as 1.084 represent exact rational numbers, not floating 
point approximations.) Up to scaling, Pf is the unique polynomial of degree 3 
such that /(O) = /(O) and pj has a double root at 22/3, and that is how it was 
constructed. 

Lemma 2.1. A contains a nonzero vector of length at most 1.084. 

Note that our lower bound (4/3)^/^ for the minimal vector length is between 
1.074 and 1.075. 

Proof. By Poisson summation, 

xga teA* 

because f{t) > for all t. On the other hand, /(O) = /(O) and f{x) < for 
\x\> 1.084, so it follows that 

E /(^) > 0- 

xeA, 0<|2:|<1.084 

Therefore the minimal vector length of A is at most 1.084. □ 
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Call a vector in A "nearly minimal" if it has length in [(4/3)^/^,1.114) (the 
reason for the upper bound of 1.114 will become clear below). 

Lemma 2.2. There are at most six nearly minimal vectors in A. 

Proof. If X and y are nearly minimal and 9 £ [0, tt] is the angle between them, then 

\^\2 ^ |^|2 _ 1^ _ y^2 2 . 1.114^ - (4/3) ^ , 

cost^ = ,, , < — , . ,„ < 0.575. 

2\x\\y\ - 2 .(4/3)1/2 

Thus, 6 > cos~i 0.575 > 2tt ■ 0.152. There is no room for seven vectors to be 
separated by such an angle: if we assign to each nearly minimal vector the arc on 
the unit circle consisting of all points within angle 27r • 0.076 of it, then these arcs 
do not overlap. Because each arc has length 27r • 0.152, if there were seven of them 
then the arc length of the unit circle would be at least 27r • 1.064. □ 

In the other direction, we can use Poisson summation to see that there are more 
than five nearly minimal vectors in A (so there must be exactly six) . First we need 
one more lemma: 

Lemma 2.3. Every nonzero vector in A is either nearly minimal or else has length 
at least 1.62. 

Proof. First, note that f{x) is a decreasing function of |a;| for |a;| G [0,1.084]. In 
the inequality 

Lemma |2 . 21 implies that there are at most six positive terms on the left hand side 
other than a; = 0. Each of them is at most /((4/3)i/'*) , where for r € [0, oo) we 
write /(r) to indicate the common value f{x) with |x| = r. Thus, 

(2.1) 6/((4/3)V4)+ J2 /(^)>0. 

xeA, f{x)<Q 

One can check that f{x) < -3/((4/3)i/'i) < for |x| e [1.114,1.62], and hence 
there is no a; g A satisfying |x| G [1.114, 1.62] (because the x and —x terms in (|2.1|) 
would combine to make the left hand side negative). □ 

Lemma 2.4. There are more than five nearly minimal vectors in A. 

Proof Define g : _^ R by g{x) = pg(27r|a;|2)e-''l^l', where 

Pg{u) = (13 - w)(1075 + 220u + 69u^). 

Then g{t) = pg(27r|t|2)e-^l*l' , where 

Pg{u) ^ (401 + 69w)(w-7)^ 

We have 5(2;) < for |x| > 1.62, and g{t) > for all t. 

It follows from Poisson summation that if A4 denotes the set of nearly minimal 
vectors in A, then 

g{0) + E gix) > g{0). 
xeM 

On the other hand, g{x) is a decreasing function of |a;| for |a;| G [(4/3)^/'^, 1.114), so 

g(0) + |M|g((4/3)i/^)>?(0). 
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Thus, 



as desired. □ 

3. Arrangement of the nearly minimal vectors 

There are two distinct nearly minimal vectors x,y £ A within angle 27r/6 of each 
other, because that is the average angle between them as one moves around the 
circle. Then x — y must be nearly minimal as well: if d is the angle between them, 
then 

|a; - yf = \xf + \yf - 2\x\\y\ cose < \xf + \yf - \x\\y\ 
because ^ < 27r/6, so 

\x-yf <2- 1.114^ - (4/3)^/^ < 1.33. 

It follows that x, y, X — y, and their negatives form the full list of nearly minimal 

vectors. 

Lemma 3.1. The vectors x and y form a basis of A- 

Proof. The nearly minimal vectors are all in the span of x and y, so we need only 
prove that the nearly minimal vectors span A. Suppose not, and that z e A is the 
smallest vector not in their span. There cannot exist a nearly minimal vector u 
such that |u — z| < \z\. For every nearly minimal u, the angle 9 between u and z 
satisfies 

cos6» = ' ' ' ,, < -LI < - = COS — 

2|m||z| - 2\z\ 2 6 

because |u — z| > \z\ and |u| < \z\. In other words, no nearly minimal vector 
is within angle 27r/6 of z. That is impossible, because then z together with the 
nearly minimal vectors would form seven vectors separated by angles of at least 
277-0.152. □ 

Thus, we have shown that A has a basis x, y such that \x\, \y\ € [(4/3)-^/*, 1.114) 
and the angle between x and y is in [2tt ■ 0.152, 2tt/Q]. 

Consider the lattice 3^/''A, which we have rescaled so that its covolume is 3^/^, 
the same as that of A2. The rescaled basis vectors 3^^^x,'S^^^y have lengths in 
[V2, 1.467) and therefore norms in [2,2.16). If 9 denotes the angle between them 
then 

(3^/'*a;,3^/''y) = 3^/^! Iz/I cos 6* G [2 cos(27r/6), 1.467^ cos(27r • 0.152)] C [1,1.243]. 

Therefore the Gram matrix of 3^/^A with respect to the basis 3^/^a;, 3^/^?/ has 
entries within 0.243 of those of the Gram matrix of A2 with respect to v, w. It is 
in this sense that A is a perturbation of A2. 

4. Local optimality 

In this section, we prove that A2 is a strict local optimum for density. It is more 
convenient to deal with quadratic forms than with lattices. Let Q be the quadratic 
form corresponding to A2 and its basis v, w (i.e., Q{s, t) = \sv + tw]"^ for s, f e M). 
The Gram matrix of A2 with respect to the basis v, w is 

'2 1^ 
,1 2> 
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and we have 




Note that the Gram matrix of a lattice basis determines the lattice up to isometrics, 
so we lose no information by focusing on it. 

The minimal norm M of Q is the minimum of Q{s,t) over {s,t) {(0,0)} 
(i.e., M = 2), and the determinant D is defined by 

D = det(j J) =3. 

In these terms, the density of A2 equals (7r/4)Z?^^/^Af, because D^l"^ is the covol- 
ume vol(M^/A2) and M is the square of the minimal vector length of A2. 

We will show that if Q is slightly perturbed, then either the perturbation is itself 
proportional to Q (which corresponds to rescaling A2) or else the density strictly 
decreases. We use symmetric perturbations because Gram matrices are symmetric. 

Let max{|a|, |c|} = p > 0, and let Qp be the perturbation of Q with matrix 

/2 + a l + 
\\ + h 2 + cj' 

Let 

Dp ^ det (^^^"^ ^tc) =3 + 2(a + c-6) + (ac-62), 
and let Mp be the minimal norm of Qp. 

Lemma 4.1. If p < 12/47, then either Dp^^^Mp < D~^/'^M or Qp is proportional 
to Q. 

Proof. We begin by assuming that a + c = b. In other words, the linear terms in 
the expansion of Dp cancel. The minimal norm of Q occurs at (1, 0), (0, 1), (1, —1), 
and their negatives, and at these points Qp takes on the values 2 + a, 2 + c, and 
2 + a — 26 + c, respectively. Consider the perturbations a, c, and a — 26 + c of these 
values away from 2. It follows from a + c = 6 that a — 26 + c = — 6 and that these 
three perturbations sum to 0. Thus, at least one of them must be negative (if all 
three vanished then p would too). In fact, one of them is at most —p/2: if a, c, or 
—6 equals —p then that is trivial, and otherwise one of them equals p and then one 
of the other two is at most —p/2 because they sum to 0. 

Thus, if a + c = 6, then Mp < 2 — p/2. When we combine that with 

Dp = i+{ac-h^) >3-2p^ 

we find that 

D-^l'^Mp < (3 - 2p^)-^/2(2 - p/2). 

When < p < 24/35, 

(3 - 2p2)-i/2(2 - pI2) < 3"^/2 • 2 = D-^^^M, 

so the density of the perturbed lattice is strictly less than that of A2 . 

In general we cannot assume a + c = b. However, Qp is proportional to a 
perturbation of Q in which that equation holds. Set A= (o— 2c+26)/(3 + a + c— 6), 
B = (46-a-c)/(3 + a + c-6), and C = (26 + c - 2a)/(3 + a + c - 6). Then 

/2 + a 1 + 6\ ^ ,..^.f2 + A 1 + B\ 

[l + b 2 + cj=(l + (« + ^-^)/^)(,l + B 2 + c) 
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and A + C = B. Clearly, \A\ < 5p/(3-3p), \B\ < 2p/(l-p), and |C| < 5p/{3-3p). 
If p < 12/47, then \A\, \B\, \C\ < 24/35, and so the perturbed lattice is either itself 
hexagonal or else strictly less dense than A2. □ 

Theorem 4.2. The hexagonal lattice A2 is the unique densest lattice in R^, up to 
scaling and isometrics o/M^. 

Proof. We have seen that if the Gram matrix of A2 is perturbed by at most 12/47, 
then either the density strictly decreases or the lattice remains the same (up to 
a similarity). On the other hand, the Gram matrix of 3"'^/^A with respect to 
3^/^y has entries within 0.243 of those of the Gram matrix of A2. Observing 
that 0.243 < 12/47 completes the proof. □ 

5. Relationship with the higher-dimensional cases 

The techniques we apply in |CK| are completely analogous to those we have 
used here, and each of the lemmas from this paper has a counterpart there. How- 
ever, there are a couple of key steps in jCKj that are not apparent from the two- 
dimensional well as some technical obstacles that must be overcome. 

5.1. Spherical codes. The biggest conceptual difference between this proof and 
the one in jCKj is in how the arrangement of nearly minimal vectors is studied. 
In each case, when rescaled to the unit sphere they form a spherical code, i.e., a 
subset of the sphere that contains no pair of points closer than some given angle. 
In the unit sphere is a circle, and distributing points on a circle with at least 
a given angle between them is trivial. In higher dimensions this problem is much 
more subtle. 

Independently, Levenshtein |Levj and Odlyzko and Sloane jOS| used linear pro- 
gramming bounds for spherical codes to solve the kissing problem in and M."^^: 
how many unit balls can be placed tangent to a given one, if they may not over- 
lap except tangentially? The answers are 240 and 196560, respectively, and the 
arrangements come from the and Leech lattice packings. 

In |CKj we use linear programming bounds in a similar way to bound the number 
of nearly minimal vectors and the angles between them. 

5.2. Spherical designs and association schemes. It is not enough in higher 
dimensions simply to compute the number of nearly minimal vectors. To deter- 
mine the configuration more precisely, we apply techniques derived from |DGS| . In 
particular, we formulate a notion of an approximate spherical design, and prove 
that the nearly minimal vectors form one. We then show that if they are grouped 
according to the (approximate) angles between them, then they form an associa- 
tion scheme. Finally, we prove that the association scheme with these parameters is 
unique, which lets us conclude that the nearly minimal vectors form a perturbation 
of the desired configuration. In principle all these results have analogues in R^, 
with similar proofs, but we did not require them for Theorem 14.21 and including 
them would have substantially lengthened this article. 

5.3. Local optimality and angle bounds. The proof of strict local optimality 
is analogous to that given here. We use Voronoi's theorem characterizing locally 
optimal lattices (as proved in |GLI §39]). However, it is more difficult to supply 
numerical bounds on how large the perturbations can be without letting the density 
increase. Section 10 of |CKj is devoted to computing such a bound. 
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In the other direction, it is also more difRcult to prove bounds on how close the 
unknown lattice is to Es or the Leech lattice. Section 7 applies the uniqueness of 
the association scheme to this problem, in a way that has no analogue in because 
fewer possible angles occur between minimal vectors in A2 than in the Leech lattice. 

5.4. Computer calculations. Completing the proofs in and R^** requires much 
sharper estimates than in M^. For example, along the way we prove that no sphere 
packing in R^"* can exceed the Leech lattice's density by a factor of 1 + 1.65 • lO"'^". 
Proving such bounds using these techniques requires applying Poisson summation 
to much more complicated functions. For example, the analogue in |CKj of the 
polynomial pf from this paper has degree 803 and rational coefficients with de- 
nominator 10'^°'^°. Checking the desired properties clearly requires a computer. We 
have arranged the calculations so as to run in one hour on an ordinary personal 
computer, and we have made our code available (see Appendix A in |('K| ). All of 
our calculations use exact rational arithmetic and are completely rigorous. 

In addition to the difficulty of checking the assertions about /, there is the issue 
of how to construct this function. We used a computer to search for polynomials 
with the desired properties, and optimized them using a high-dimensional version of 
Newton's method. The whole process took our computers approximately a month, 
but fortunately one can check the proof simply by verifying the final answer, with 
no need to reconstruct it from scratch. 

6. Future prospects 

There is no likelihood that our techniques will work in any dimension except 1, 
2, 8, or 24. They depend on using the methods of jCE| to prove nearly sharp upper 
bounds for the packing density, and that does not seem to happen in any other 
dimension. It seems counterintuitive that determining the densest lattice in ]S}^ 
appears much more difficult than in R^ or R^"*, but that simply seems to be the 
case. 

One might hope to solve the full sphere packing problem using our techniques. 
We cannot absolutely rule that out, but we do not consider it feasible. There are 
at least two substantial obstacles: 

(1) It is not known whether or the Leech lattice is even locally optimal when 
the perturbations leave the space of lattices. 

(2) In a general sphere packing, different spheres can be tangent to different 
numbers of other spheres (unlike the case of lattices). Our method of prov- 
ing lower bounds for the number of nearly minimal vectors can be used to 
bound the average number of near neighbors of a sphere in a dense packing. 
However, we cannot rule out the possibility of a small minority of spheres 
with few neighbors. 

We do not know how to deal with the first problem, but perhaps it could be 
done. The second problem appears quite fundamental, and we suspect that dealing 
with it would require a major advance in the theory of sphere packing. 

Instead, we are convinced that the right way to solve the sphere packing problem 
in R^ and R^"* is to prove the following conjecture from |CE| . Call a continuous 
function / : R" R admissible if there is a constant <5 > such that both |/(a;)| 
and \f{x)\ are bounded above by a constant times (1 + |a;|)~"~''. (Note that the 
bound on / implies that the integral defining / converges.) 
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Conjecture 6.1. For {n,r) = (8,\/2) or (24,2), there is an admissible function 
/ : R" -> M such that /(O) = /(O) = 1, f{x) < for \x\ > r, and J{t) > for all t. 

If Coniecture 16 . II holds . then Theorem 3.2 of |UE| implies that Es and the Leech 
lattice are the densest sphere packings in and K^'', respectively. In |CK| we 
achieve r < 2(l + 6.851 • 10^^^) when n = 24, which is numerical evidence that 
r = 2 could be achieved. 

Acknowledgements 
We thank Amanda Beeson and James Bernhard for comments on our paper. 

Appendix A. Poisson summation 

To make this paper self-contained, we provide here the standard proof of the 
Poisson summation formula: 

Theorem A.l (Poisson summation). Let f : M" ^ R 6e a Schwartz function and 
A C R" a lattice. Then 

^ ^("'^ ^ vol(RVA) ^ 
Proof. It is easier to prove a more general formula. Define F : R" — > R by 

F(z) = ^/(a: + z). 

a:6A 

In other words, F is f made periodic modulo A, and we wish to compute F{Q). 
Because / is a Schwartz function, the sum defining F converges and defines a C°° 
function. It follows that we can expand F as a Fourier series, in particular as a 
linear combination of the exponential functions that are periodic modulo A. The 
function 

is periodic modulo A iff t G A*. Thus, 

F{z) = J2 cte-'"*<*'^> 
teA* 

for some coefficients q. We can compute Ct using orthogonality: if I? is a funda- 
mental domain for A, then 



ID 

When we substitute vol(Z3) = vol(R"/A) and the definition of F{z) into this equa- 
tion, we find that 



Ct 



,^ A JD + x 

/(z)e2"<*'^) dz 



vol( 


;r"/a) 




1 


vol( 


;r"/a) 




1 


vol( 


;r"/a) 




1 



vol(R"/A) 
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Here D + x = {d + x : d £ D}, and these sets tile R" as x ranges over A. 
It follows that 

xeA ^ ' ' teA* 

Setting z = yields the desired result. □ 
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